Insight into small-molecule inhibitors targeting extracellular nucleotide pyrophosphatase/phosphodiesterase1 for potential multiple human diseases.

Extracellular nucleotide pyrophosphatase/phosphodiesterase 1 (ENPP1) has been identified as a type II transmembrane glycoprotein. It plays a crucial role in various biological processes, such as bone mineralization, cancer cell proliferation, and immune regulation. Consequently, ENPP1 has garnered attention as a promising target for pharmacological interventions. Despite its potential, the development of clinical-stage ENPP1 inhibitors for solid tumors, diabetes, and silent rickets remains limited. However, there are encouraging findings from preclinical trials involving small molecules exhibiting favorable therapeutic effects and safety profiles. This perspective aims to shed light on the structural properties, biological functions and the relationship between ENPP1 and diseases. Additionally, it focuses on the structure-activity relationship of ENPP1 inhibitors, with the intention of guiding the future development of new and effective ENPP1 inhibitors.

Differential effects of the lipidic and ionic microenvironment on NPP1's phosphohydrolase and phosphodiesterase activities.

Ecto-nucleotide pyrophosphatase/phosphodiesterase 1 (NPP1) is an enzyme present in matrix vesicles (MV). NPP1 participates on the regulation of bone formation by producing pyrophosphate (PPi) from adenosine triphosphate (ATP). Here, we have used liposomes bearing dipalmitoylphosphatidylcholine (DPPC), sphingomyelin (SM), and cholesterol (Chol) harboring NPP1 to mimic the composition of MV lipid rafts to investigate ionic and lipidic influence on NPP1 activity and mineral propagation. Atomic force microscopy (AFM) revealed that DPPC-liposomes had spherical and smooth surface. The presence of SM and Chol elicited rough and smooth surface, respectively. NPP1 insertion produced protrusions in all the liposome surface. Maximum phosphodiesterase activity emerged at 0.082 M ionic strength, whereas maximum phosphomonohydrolase activity arose at low ionic strength. Phosphoserine-Calcium Phosphate Complex (PS-CPLX) and amorphous calcium-phosphate (ACP) induced mineral propagation in DPPC- and DPPC:SM-liposomes and in DPPC:Chol-liposomes, respectively. Mineral characterization revealed the presence of bands assigned to HAp in the mineral propagated by NPP1 harbored in DPPC-liposomes without nucleators or in DPPC:Chol-liposomes with ACP nucleators. These data show that studying how the ionic and lipidic environment affects NPP1 properties is important, especially for HAp obtained under controlled conditions in vitro.

A Brief Overview of the Toxic Sphingomyelinase Ds of Brown Recluse Spider Venom and Other Organisms and Simple Methods To Detect Production of Its Signature Cyclic Ceramide Phosphate.

A special category of phospholipase D (PLD) in the venom of the brown recluse spider (Loxosceles reclusa) and several other sicariid spiders accounts for the dermonecrosis and many of the other clinical symptoms of envenomation. Related proteins are produced by other organisms, including fungi and bacteria. These PLDs are often referred to as sphingomyelinase Ds (SMase Ds) because they cleave sphingomyelin (SM) to choline and "ceramide phosphate." The lipid product has actually been found to be a novel sphingolipid: ceramide 1,3-cyclic phosphate (Cer1,3P). Since there are no effective treatments for the injury induced by the bites of these spiders, SMase D/PLDs are attractive targets for therapeutic intervention, and some of their features will be described in this minireview. In addition, two simple methods are described for detecting the characteristic SMase D activity using a fluorescent SM analog, (N-[12-[(7-nitro-2-1,3-benzoxadiazol-4-yl)amino]dodecanoyl]-SM (C12-NBD-SM), that is cleaved to C12-NBD-Cer1,3P, which is easily separated from other potential metabolites by thin-layer chromatography and visualized under UV light. Besides confirming that C12-NBD-Cer1,3P is the only product detected upon incubation of C12-NBD-SM with brown recluse spider venom, the method was also able to detect for the first time very low levels of activity in venom from another spider, Kukulcania hibernalis The simplicity of the methods makes it relatively easy to determine this signature activity of SMase D/PLD. SIGNIFICANCE STATEMENT: The sphingomyelinase D/phospholipase D that are present in the venom of the brown recluse spider and other sources cause considerable human injury, but detection of the novel sphingolipid product, ceramide 1,3-cyclic phosphate, is not easy by previously published methods. This minireview describes simple methods for detection of this activity that will be useful for studies of its occurrence in spider venoms and other biological samples, perhaps including lesions from suspected spider bites and infections.

The combined treatment of gold nanoparticles associated with photobiomodulation accelerate the healing of dermonecrotic lesion.

Introduction: The search for fast and efficient treatment for dermonecrotic lesions caused by the venom of the spider from the Loxosceles simillis, is a demand in health. Prednisolone is one of the most used drugs, however it has side effects. In this context, addictionally gold nanoparticles (GNPs) have anti-inflammatory, antioxidant, and antibacterial properties. The use of photobiomodulation has show to be efficient in the process of tissue repair. Therefore, the purpose of this study was to investigate the anti-inflammatory effect of photobiomodulation and GNPs associated or not with a low concentration of prednisolone in animal models of dermonecrotic lesion.Methodology: For this, rabbits with venon-induced dermonecrotic lesion were subjected to topical treatment with prednisolone + laser or GNPs + laser or Pred-GNPs + laser. The area of edema, necrosis and erythema were measured. On the last day of treatment, the animals were euthanized to remove the organs for histopathological and biochemical analysis.Results: All treatments combinations were effective in promoting the reduction of necrotic tissue and erythema.Conclusion: With this results, we suggest that the use of laser and nanoparticles, associated or not with prednisolone, should be considered for the treatment of dermonecrotic injury.

Synthesis and degradation of the cyclic dinucleotide messenger c-di-AMP in the hyperthermophilic archaeon Pyrococcus yayanosii.

Cyclic di-adenosine monophosphate (c-di-AMP) is a newly identified prokaryotic cyclic dinucleotide second messenger well elucidated in bacteria, while less studied in archaea. Here, we describe the enzymes involved in c-di-AMP metabolism in the hyperthermophilic archaeon Pyrococcus yayanosii. Our results demonstrate that c-di-AMP is synthesized from two molecules of ATP by diadenylate cyclase (DAC) and degraded into pApA and then to AMP by a DHH family phosphodiesterase (PDE). DAC can be activated by a wider variety of ions, using two conserved residues, D188 and E244, to coordinate divalent metal ions, which is different from bacterial CdaA and DisA. PDE possesses a broad substrate spectrum like bacterial DHH family PDEs but shows a stricter base selection between A and G in cyclic dinucleotides hydrolysis. PDE shows differences in substrate binding patches from bacterial counterparts. C-di-AMP was confirmed to exist in Thermococcus kodakarensis cells, and the deletion of the dac or pde gene supports that the synthesis and degradation of c-di-AMP are catalyzed by DAC and PDE, respectively. Our results provide a further understanding of the metabolism of c-di-AMP in archaea.

Isolation and characterization of the first phosphodiesterase (Bj-PDE) from the venom of Bothrops jararacussu snake.

Phosphodiesterases are exonucleases that sequentially hydrolyse phosphodiester bonds of polynucleotides from the 3'-end and release 5-mononucleotides. After more than one decade without any advance in the study of Bothropic phosphodiesterases, we described here the isolation of the first phosphodiesterase from Bothrops jararacussu, which we named Bj-PDE. A five-step column chromatography procedure (size exclusion, hydrophobic interaction, cation exchange, lentil lectin affinity, and blue sepharose affinity) enabled isolation of Bj-PDE with preserved and stable enzymatic activity (bis(p-nitrophenyl) phosphate substrate), Km = 6.9 mM (± 0.7 mM), kcat/Km = 1.7 × 104 M-1 s-1 (± 0.2 × 104 M-1 s-1), MW = 116 kDa (SDS-PAGE), optimum activity around 45 °C at pH 8.0, and stability for 81 days at different storage temperatures (8, -20, and - 80 °C). Ca2+ and Mg2+ ions positively influenced Bj-PDE activity, while EDTA had the opposite action. Zn2+ restored >50 % of enzyme activity after its inhibition by EDTA. The Bj-PDE partial sequence identified by mass spectrometry was very similar to the sequence of BATXPDE1 from Bothrops atrox, which was evolutionarily close to this new PDE. Therefore, our study represents an important progress on the isolation of this minor toxin and sheds new lights on the properties and bioprospection of bothropic phosphodiesterases.

Brown Spider Venom Phospholipase-D Activity upon Different Lipid Substrates.

Brown spider envenomation results in dermonecrosis, characterized by an intense inflammatory reaction. The principal toxins of brown spider venoms are phospholipase-D isoforms, which interact with different cellular membrane components, degrade phospholipids, and generate bioactive mediators leading to harmful effects. The Loxosceles intermedia phospholipase D, LiRecDT1, possesses a loop that modulates the accessibility to the active site and plays a crucial role in substrate. In vitro and in silico analyses were performed to determine aspects of this enzyme's substrate preference. Sphingomyelin d18:1/6:0 was the preferred substrate of LiRecDT1 compared to other Sphingomyelins. Lysophosphatidylcholine 16:0/0:0 was preferred among other lysophosphatidylcholines, but much less than Sphingomyelin d18:1/6:0. In contrast, phosphatidylcholine d18:1/16:0 was not cleaved. Thus, the number of carbon atoms in the substrate plays a vital role in determining the optimal activity of this phospholipase-D. The presence of an amide group at C2 plays a key role in recognition and activity. In silico analyses indicated that a subsite containing the aromatic residues Y228 and W230 appears essential for choline recognition by cation-π interactions. These findings may help to explain why different cells, with different phospholipid fatty acid compositions exhibit distinct susceptibilities to brown spider venoms.

Autotaxin and Lysophosphatidic Acid Signalling: the Pleiotropic Regulatory Network in Cancer.

Autotaxin, also known as ecto-nucleotide pyrophosphatase/phosphodiesterase family member 2, is a secreted glycoprotein that plays multiple roles in human physiology and cancer pathology. This protein, by converting lysophosphatidylcholine into lysophosphatidic acid, initiates a complex signalling cascade with significant biological implications. The article outlines the autotaxin gene and protein structure, expression regulation and physiological functions, but focuses mainly on the role of autotaxin in cancer development and progression. Autotaxin and lysophosphatidic acid signalling influence several aspects of cancer, including cell proliferation, migration, metastasis, therapy resistance, and interactions with the immune system. The potential of autotaxin as a diagnostic biomarker and promising drug target is also examined.

From Myricetin to the Discovery of Novel Natural Human ENPP1 Inhibitors: A Virtual Screening, Molecular Docking, Molecular Dynamics Simulation, and MM/GBSA Study.

It was recently revealed that naturally occurring myricetin can inhibit ectonucleotidase ectonucleotide pyrophosphatase/phosphodiesterase 1 (ENPP1), which, in turn, can treat ischemic cardiac injury. However, due to myricetin's poor druggability, its further developments are relatively limited, which necessitates the discovery of novel ENPP1-inhibiting myricetin analogs as alternatives. In this study, the binding model of myricetin with ENPP1 was elucidated by molecular docking and molecular dynamics studies. Subsequently, virtual screening on the self-developed flavonoid natural product database (FNPD), led to the identification of two flavonoid glycosides (Cas No: 1397173-50-0 and 1169835-58-8), as potential ENPP1 inhibitors. Docking scores and MM/GBSA binding energies predicted that they might have higher inhibitory effects than myricetin. This study provides a strong foundation for the future development of ischemic cardiac injury drugs.

Twin cyclic mononucleotide cyclase and phosphodiesterase domain architecture as a common feature in complex plant proteins.

The majority of proteins in both prokaryote and eukaryote proteomes consist of two or more functional centers, which allows for intramolecular tuning of protein functions. Such architecture, as opposed to animal orthologs, applies to the plant cyclases (CNC) and phosphodiesterases (PDEs), the vast majority of which are part of larger multifunctional proteins. In plants, until recently, only two cases of combinations of CNC-PDE in one protein were reported. Here we propose that in plants, multifunctional proteins in which the PDE motif has been identified, the presence of the additional CNC center is common. Searching the Arabidopsis thaliana proteome with a combined PDE-CNC motif allowed the creation of a database of proteins with both activities. One such example is methylenetetrahydrofolate dehydrogenase, in which we determined the activities of adenylate cyclase (AC) and PDE. Based on biochemical and mutagenesis analyses we assessed the impact of the AC and PDE catalytic centers on the dehydrogenase activity. This allowed us to propose additional regulatory mechanism that govern folate metabolism by cAMP. It is therefore conceivable that the combined CNC-PDE architecture is a common regulatory configuration, where control of the level of cyclic nucleotides (cNMP) influences other catalytic activities of the protein.

Artificial phosphatase upon premicellar nanoarchitectonics of lanthanum complexes with long-chained imidazole derivatives.

Four novel long chain-containing tridentate imidazole derivatives (Ln, n = 1, 2, 3, 4) were synthesized for in situ formation of mononuclear lanthanum(III) complexes as artificial phosphodiesterases. These in-situ formed La(III) complexes (named LaLn) were used to catalyze the transesterification of 2-hydroxypropyl p-nitrophenyl phosphate (HPNP), a classic RNA model. Critical aggregation concentrations (CAC) were determined for the as-prepared tridentate imidazole derivatives as ligands and corresponding mixtures of equivalent ligand and La3+ ion with a mole rate of 1:1. It denotes that the introduction of La3+ ion increases the CAC values of imidazole derivatives by about 2 to 3 folds. Foaming test shows that the foam height is positively correlated with the length of hydrophobic chain. Transesterification of HPNP mediated by LaLn nanoarchitectonics indicates that the introducing of hydrophobic chain benefits rate enhancement, showing excess three orders of magnitude acceleration under physiological conditions (pH 7.0, 25 °C). Moreover, catalytic reactivities of these La(III) complexes increased along with the increase in chain length: LaL1 < LaL2 < LaL3 < LaL4, suggesting a positive correlation to hydrophobic chain length.

Vibrio cholerae V-cGAP3 Is an HD-GYP Phosphodiesterase with a Metal Tunable Substrate Selectivity.

Cyclic dinucleotides (CDNs) are signaling molecules involved in the immune response and virulence factor production. CDN cellular levels are fine-tuned by metal-dependent phosphodiesterases (PDEs), among which HD-GYPs make up a subclass of the larger HD-domain protein superfamily. The human pathogen Vibrio cholerae (Vc) encodes nine HD-GYPs, one of which is V-cGAP3 (or VCA0931). V-cGAP3 acts on c-di-GMP and 3'3'c-GAMP, and this activity is related to bacterial infectivity. However, the extant chemical makeup of the V-cGAP3 cofactor and steady state parameters have not been established. Employing electron paramagnetic resonance and Mössbauer spectroscopy in tandem with elemental analyses and activity assays, we demonstrate that V-cGAP3 coordinates different dimetal cofactors with variable activities. MnII and FeII afford c-di-GMP hydrolysis with the highest observed rates, while c-GAMP hydrolysis is selectively dependent on Mn. V-cGAP3 has a single functional domain, and this simple architecture allows us to examine the roles of characteristic conserved residues in catalysis. Substitution of the adjacent to the active site GYP residue triad and the specifically conserved in HD-domain PDEs fifth histidine ligand (i.e., H371 in V-cGAP3) with alanines severely compromises CDN hydrolysis but only modestly affects cofactor incorporation. Our data are consistent with V-cGAP3 being the major regulator of 3'3'c-GAMP hydrolysis in Vc and delineate the importance of specific residues in tuning activity in HD-GYPs in general. We propose that HD-GYPs exhibit diversity in their metallocofactors and substrates, which may serve to increase their functional potential in regulatory pathways or allow for PDE activity upon adaptation of the parent organism to diverse environmental niches.

In Vitro and In Silico Studies of Human Tyrosyl-DNA Phosphodiesterase 1 (Tdp1) Inhibition by Stereoisomeric Forms of Lipophilic Nucleosides: The Role of Carbohydrate Stereochemistry in Ligand-Enzyme Interactions.

Inhibition of human DNA repair enzyme tyrosyl-DNA phosphodiesterase 1 (Tdp1) by different chiral lipophilic nucleoside derivatives was studied. New Tdp1 inhibitors were found in the series of the studied compounds with IC50 = 2.7-6.7 µM. It was shown that D-lipophilic nucleoside derivatives manifested higher inhibition activity than their L-analogs, and configuration of the carbohydrate moiety can influence the mechanism of Tdp1 inhibition.

Structural and functional identification of the uncharacterized metallo-ß-lactamase superfamily protein TW9814 as a phosphodiesterase with unique metal coordination.

Metallo-ß-lactamase (MBL) superfamily proteins have a common αß/ßα sandwich fold and perform a variety of functions through metal-mediated catalysis. However, because of the enormous scale of this superfamily, only a small percentage of the proteins belonging to the superfamily have been annotated structurally or functionally to date. Therefore, much remains unknown about the MBL superfamily proteins. Here, TW9814, a hypothetical MBL superfamily protein, was structurally and functionally investigated. Guided by the crystal structure of dimeric TW9814, it was demonstrated that TW9814 functions as a phosphodiesterase (PDE) in the presence of divalent metal ions such as manganese(II) or nickel(II). A docking model between TW9814 and the substrate bis(p-nitrophenyl)phosphate (bpNPP) showed the importance of the dimerization of TW9814 for its bpNPP-hydrolyzing activity and for the interaction between the enzyme and the substrate. TW9814 showed outstanding catalytic efficiency (kcat/Km) under alkaline conditions compared with other PDEs. The activity of TW9814 appears to be regulated through a disulfide bond, which is a feature that is not present in other MBL superfamily members. This study provides a platform for the functional characterization of other hypothetical proteins of the MBL or other superfamilies.

Pyranodipyran Derivatives with Tyrosyl DNA Phosphodiesterase 1 Inhibitory Activities and Fluorescent Properties from Aspergillus sp. EGF 15-0-3.

Four new benzodipyran racemates, namely (±)-aspergiletals A-D (3-6), representing a rare pyrano[4,3-h]chromene scaffold were isolated together with eurotiumide G (1) and eurotiumide F (2) from the soft-coral-derived fungus Aspergillus sp. EGF 15-0-3. All the corresponding optically pure enantiomers were successfully separated by a chiral HPLC column. The structures and configurations of all the compounds were elucidated based on the combination of NMR and HRESIMS data, chiral separation, single-crystal X-ray diffraction, quantum chemical 13C NMR, and electronic circular dichroism calculations. Meanwhile, the structure of eurotiumide G was also revised. The TDP1 inhibitor activities and photophysical properties of the obtained compounds were evaluated. In the TDP1 inhibition assay, as a result of synergy between (+)-6 and (-)-6, (±)-6 displayed strong inhibitory activity to TDP1 with IC50 values of 6.50 ± 0.73 µM. All compounds had a large Stokes shift and could be utilized for elucidating the mode of bioactivities by fluorescence imaging.

Specificity of Loxosceles α clade phospholipase D enzymes for choline-containing lipids: Role of a conserved aromatic cage.

Spider venom GDPD-like phospholipases D (SicTox) have been identified to be one of the major toxins in recluse spider venom. They are divided into two major clades: the α clade and the ß clade. Most α clade toxins present high activity against lipids with choline head groups such as sphingomyelin, while activities in ß clade toxins vary and include preference for substrates containing ethanolamine headgroups (Sicarius terrosus, St_ßIB1). A structural comparison of available structures of phospholipases D (PLDs) reveals a conserved aromatic cage in the α clade. To test the potential influence of the aromatic cage on membrane-lipid specificity we performed molecular dynamics (MD) simulations of the binding of several PLDs onto lipid bilayers containing choline headgroups; two SicTox from the α clade, Loxosceles intermedia αIA1 (Li_αIA) and Loxosceles laeta αIII1 (Ll_αIII1), and one from the ß clade, St_ßIB1. The simulation results reveal that the aromatic cage captures a choline-headgroup and suggest that the cage plays a major role in lipid specificity. We also simulated an engineered St_ßIB1, where we introduced the aromatic cage, and this led to binding with choline-containing lipids. Moreover, a multiple sequence alignment revealed the conservation of the aromatic cage among the α clade PLDs. Here, we confirmed that the i-face of α and ß clade PLDs is involved in their binding to choline and ethanolamine-containing bilayers, respectively. Furthermore, our results suggest a major role in choline lipid recognition of the aromatic cage of the α clade PLDs. The MD simulation results are supported by in vitro liposome binding assay experiments.

Structure-based discovery of (S)-2-amino-6-(4-fluorobenzyl)-5,6,11,11a-tetrahydro-1H-imidazo[1',5':1,6]pyrido[3,4-b]indole-1,3(2H)-dione as low nanomolar, orally bioavailable autotaxin inhibitor.

Inhibition of extracellular secreted enzyme autotaxin (ATX) represents an attractive strategy for the development of new therapeutics to treat various diseases and a few inhibitors entered in clinical trials. We herein describe structure-based design, synthesis, and biological investigations revealing a potent and orally bioavailable ATX inhibitor 1. During the molecular docking and scoring studies within the ATX enzyme (PDB-ID: 4ZGA), the S-enantiomer (Gscore = -13.168 kcal/mol) of the bound ligand PAT-494 scored better than its R-enantiomer (Gscore = -9.562 kcal/mol) which corroborated with the reported observation and analysis of the results suggested the scope of manipulation of the hydantoin substructure in PAT-494. Accordingly, the docking-based screening of a focused library of 10 compounds resulted in compound 1 as a better candidate for pharmacological studies. Compound 1 was synthesized from L-tryptophan and evaluated against ATX enzymatic activities with an IC50 of 7.6 and 24.6 nM in biochemical and functional assays, respectively. Further, ADME-PK studies divulged compound 1 as non-cytotoxic (19.02% cell growth inhibition at 20 µM in human embryonic kidney cells), metabolically stable against human liver microsomes (CLint  = 15.6 µl/min/mg; T1/2  = 113.2 min) with solubility of 4.82 µM and orally bioavailable, demonstrating its potential to be used for in vivo experiments.

Structure and function of the ecto-nucleotide pyrophosphatase/phosphodiesterase (ENPP) family: Tidying up diversity.

Ecto-nucleotide pyrophosphatase/phosphodiesterase (ENPP) family members (ENPP1-7) have been implicated in key biological and pathophysiological processes, including nucleotide and phospholipid signaling, bone mineralization, fibrotic diseases, and tumor-associated immune cell infiltration. ENPPs are single-pass transmembrane ecto-enzymes, with notable exceptions of ENPP2 (Autotaxin) and ENNP6, which are secreted and glycosylphosphatidylinositol (GPI)-anchored, respectively. ENNP1 and ENNP2 are the best characterized and functionally the most interesting members. Here, we review the structural features of ENPP1-7 to understand how they evolved to accommodate specific substrates and mediate different biological activities. ENPPs are defined by a conserved phosphodiesterase (PDE) domain. In ENPP1-3, the PDE domain is flanked by two N-terminal somatomedin B-like domains and a C-terminal inactive nuclease domain that confers structural stability, whereas ENPP4-7 only possess the PDE domain. Structural differences in the substrate-binding site endow each protein with unique characteristics. Thus, ENPP1, ENPP3, ENPP4, and ENPP5 hydrolyze nucleotides, whereas ENPP2, ENPP6, and ENNP7 evolved as phospholipases through adaptions in the catalytic domain. These adaptations explain the different biological and pathophysiological functions of individual members. Understanding the ENPP members as a whole advances our insights into common mechanisms, highlights their functional diversity, and helps to explore new biological roles.

Characterization of recombinant murine GDE4 and GDE7, enzymes producing lysophosphatidic acid and/or cyclic phosphatidic acid.

GDE4 and GDE7 are membrane-bound enzymes that exhibit lysophospholipase D activities. We found that GDE7 produced not only lysophosphatidic acid (LPA) but also cyclic phosphatidic acid (cPA) from lysophospholipids by a transphosphatidylation reaction. In contrast, GDE4 produced only LPA. The analysis of substrate specificity showed that 1-alkyl-lysophosphospholipids were preferred substrates for both enzymes rather than 1-alkyl-lysophospholipids and 1-alkenyl-lysophospholipids. Among the various lysophospholipids with different polar head groups that were tested, lysophosphatidylglycerol and lysophosphatidylserine were preferred substrates for GDE4 and GDE7, respectively. The detailed analysis of the dependency of the enzyme activities of GDE4 and GDE7 on divalent cations suggested multiple divalent cations were bound in the active sites of both enzymes. Taken together, these results suggest the possibility that GDE7 functions as a cPA-producing enzyme in the body.

Selectivity mechanism of phosphodiesterase isoform inhibitor through in silico investigations.

Understanding the selectivity mechanism of inhibitors towards homology proteins helps to design selective candidates. Phosphodiesterase (PDE) family members act in the degradation of cAMP and cGMP, among which some isoforms such as PDE9A are attracting interest for Alzheimer's disease treatment, while PDE10A is used as target for treating schizophrenia. In this study, computational methods were used to investigate the major features of PDE9A/10A, with the purpose to provide deep understanding of the molecular mechanism of selective inhibition towards these two isoforms. Our result revealed that two conserved residues Gln453 and Phe456 were proven to be crucial for the binding affinity and inhibitory selectivity of PDE9A inhibitors. In addition, the high-affinity PDE9A inhibitors always interact with the conservative hydrophobic pocket as well as Tyr424 and Ala452 of PDE9A, while PDE10A selective inhibitors need to have two hydrophobic groups and two hydrogen bond donors to interact with the conservative Tyr693, Gln726, and Phe729 of PDE10A. This study provides valuable insights into the underlying mechanism of selective inhibition targeting PDE9A and PDE10A, for further search for potent and highly selective PDE9A/10A inhibitors.